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Abstract 

Ca 2+ release into the cytosol through inositol 1,4,5-trisphosphate receptors (IP 3 Rs) plays a relevant role in numerous 
physiological processes. IP 3 R-mediated Ca 2+ signals involve Ca 2+ -induced Ca 2+ -release (CICR) whereby Ca 2+ release 
through one open IP 3 R induces the opening of other channels. IP 3 Rs are apparently organized in clusters. The signals can 
remain localized (i.e., Ca 2+ puffs) if CICR is limited to one cluster or become waves that propagate between clusters. Ca 2+ 
puffs are the building blocks of Ca 2+ waves. Thus, there is great interest in determining puff properties, especially in view of 
the current controversy on the spatial distribution of activatable IP 3 Rs. Ca 2+ puffs have been observed in intact cells with 
optical techniques proving that they are intrinsically stochastic. Obtaining a correct picture of their dynamics then entails 
being able to detect the whole range of puff sizes. Ca 2+ puffs are observed using visible single-wavelength Ca 2+ dyes, 
slow exogenous buffers (e.g., EGTA) to disrupt inter-cluster CICR and UV-photolyzable caged IP 3 . Single-wavelength dyes 
increase their fluorescence upon calcium binding producing images that are strongly dependent on their kinetic, transport 
and photophysical properties. Determining the artifacts that the imaging setting introduces is particularly relevant when 
trying to analyze the smallest Ca 2+ signals. In this paper we introduce a method to estimate the expected signal-to-noise 
ratio of Ca 2+ imaging experiments that use single-wavelength dyes. The method is based on the Number and Brightness 
technique. It involves the performance of a series of experiments and their subsequent analysis in terms of a fluorescence 
fluctuation model with which the model parameters are quantified. Using the model, the expected signal-to-noise ratio is 
then computed. Equivalence classes between different experimental conditions that produce images with similar signal-to- 
noise ratios can then be established. The method may also be used to estimate the smallest signals that can reliably be 
observed with each setting. 
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Introduction 

Calcium signals are ubiquitous [1]. Their versatility relies on the 
variety of spatio-temporal behaviors that the intracellular calcium 
concentration can display. Ca 2+ release into the cytosol through 
inositol 1,4,5-trisphosphate receptors (IP 3 Rs) is a key component 
of the Ca 2+ signaling toolkit [2]. Optical techniques and Ca 2+ 
fluorescent dyes provide a relatively non-invasive means to study 
the dynamics of intracellular Ca 2+ signals [3,4], in particular, 
those that are IP 3 R-mediated. These observations have revealed a 
wide variety of signals that go from those that remain spatially 
localized to those that propagate throughout the cell [5]. The 
observations, however, are indirect. Most Ca 2+ dyes change their 
spectral properties upon Ca 2+ binding. Thus, the observed 
fluorescence is related to the Ca 2 + -bound dye rather than the 
free Ca 2 + concentration. Having reliable estimates of the dye 
physical properties is then necessary to quantitate the underlying 
free Ca 2 + distribution [6]. IP 3 R-mediated Ca 2+ signals are 
observed using single-wavelength dyes that are excited with visible 



light and caged IP3 that is photolyzed with UV illumination to 
evoke the signals [7]. Differendy from ratiometric dyes, single- 
wavelength dyes do not allow for a direct measurement of the 
Ca 2+ concentration [8,9]. Fluorescence variations during 
Ca + signals are then presented as ratios with respect to basal 
fluorescence (prior to signal evocation) when using single- 
wavelength dyes. Namely, they are shown in terms of the ratio: 

i r u F( Xi ,tj)-Fo(xi) 

where F(Xi,tj) is the fluorescence at the pixel identified by the 
position Xj and time tj and Fq(xi) is the fluorescence at spatial 
point, Xi, averaged over time prior to signal evocation. This 
minimizes the artifacts that spatial heterogeneities due to uneven 
dye distribution, specimen thickness or illumination intensity can 
introduce and allows a direct comparison across the image. The 
question naturally arises of whether this is enough to compare 
images performed under different experimental conditions. More 
specifically, how similarly a given underlying free 
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Ca + distribution is reflected in images that are obtained with 
different dyes and/or dye concentrations or with different 
experimental set-ups. This question is particularly relevant when 
the underlying dynamics is subject to Calcium Induced Calcium 
Release (CICR) {i.e., when the Ca 2 + released through one open 
channel induces the opening of neighboring ones) since the dye or 
other substances introduced in the cells during the experiments 
can interfere with signal propagation. 

IP :i R-mediated Ca 2+ signals are subject to CICR. Namely, 
IP :i Rs need to bind IP3 and Ca 2 + to become open [10]. In most 
cells IP 3 Rs are organized in clusters. Thus IP 3 R-mediated 
Ca 2 + signals remain localized (e.g., puffs that involve the release 
through a few IP 3 Rs in a cluster) or global (e.g., waves) depending 
on CICR efficiency. Puffs are the building blocks of global signals 
and, as such, have been the subject of numerous studies [1 1]. Puffs 
are highly stochastic. On one hand, not all clusters contain the 
same number or spatial distribution of IP3RS [3,12]. On the other 
hand, the same cluster can give rise to puffs of different sizes 
depending on the number and spatial distribution of the 
activatable IP :i Rs that there are at the beginning of the signal 
[13], two properties that change with time as IP 3 binds/unbinds 
IP 3 Rs and IP3RS enter or leave their inhibited state [14]. In order 
to obtain a good understanding of the dynamics of puffs it is then 
necessary to collect enough information so as to derive an accurate 
statistical description. The latter depends on the ability of the 
experimental set-up to detect most of the evoked events. The 
smallest ones, however, can easily go undetected depending on 
how large the change in fluorescence they produce is compared to 
spontaneous fluorescence fluctuations. Being able to detect the 
smallest possible events is particularly relevant in the case of IP 3 Rs 
since there is some controversy on what their actual spatial 
distribution is and how it can change with IP 3 stimulation [15,16]. 
In order to collect enough statistics on puffs, experiments are 
performed introducing a slow exogenous Ca 2+ buffer (typically, 
EGTA) in the cell. In this way the evoked signals can be forced to 
remain local. The problem is that EGTA not only changes the 
dynamics of the signals. In principle, it can also modify the way 
that an underlying Ca 2+ distribution is reported by the observed 
fluorescence. In this paper we introduce a method to estimate the 
signal-to-noise ratio of Ca 2+ imaging experiments that use single- 
wavelength Ca 2+ dyes and EGTA. The method relies on a model 
of fluctuations that can be quantitated by means of simple 
experiments. Once the model is quantified, numerically simulated 
Ca 2 + images with fluctuations that are in accordance with the used 
experimental setting can be generated. The quantified model can 
then be used to probe the limits of signal detectability of the 
experiments. 

The method we introduce in this paper is based on the Number 
and Brightness (N&B) technique of [17]. N&B serves to separate 
the contributions of two sources of fluctuations that affect 
fluorescence images: variations in the number of molecules that 
contribute to the fluorescence at each pixel and variations in the 
number of detected photons. The N&B method has been applied 
mainly to analyze static images of samples containing fluorescent 
proteins and identify whether increments in the observed 
fluorescence are due to a change in the number of emitting 
molecules or in the brightness (i.e., the number of detected 
photons). To achieve this goal, images of the marked sample are 
obtained for different illumination powers. The mean and the 
variance of the fluorescence are then computed for each image. 
The analysis of how these quantities change with the illumination 
power allows a separation between number and brightness and 
between the mobile and the immobile fraction of the fluorescent 



molecules. In this way a pixel by pixel map of molecular number 
and aggregation can be obtained for the image. The aim of our 
method is to have a quantitative description of the distribution of 
fluorescence fluctuations when the emitters are single-wavelength 
Ca 2+ dyes in the presence of Ca 2+ . In this case, under stationary 
conditions (i.e., in the absence of Ca 2+ signals), there is an 
additional source of fluctuations with respect to those considered 
in the N&B method. Namely, single wavelength Ca 2+ dyes emit in 
the same range of wavelengths if the dye molecules are bound or 
not toCa 2+ , albeit, with a different quantum efficiency. Ca 2+ 
binding to the dye is a dynamic process. Thus, images obtained in 
the cytosol of cells containing one such dye, under stationary 
conditions, are also affected by variations in the fraction of Ca 2+ - 
bound dye molecules that contribute to the fluorescence at each 
pixel. Furthermore, this fraction depends on the number of free 
Ca 2+ ions in the region of interest which is also a dynamic 
variable. Thus, the number of quantifiable parameters of the 
fluctuation model is larger in our method than in N&B and then 
requires more experiments. The purpose of our method, however, 
is more restricted than that of the N&B technique. We are not 
trying to establish a map of the number of dye molecules of Ca 2+ 
ions that are in the region that contribute to the fluorescence at 
each pixel. We are just interested in having a quantified version of 
the distribution of fluorescence fluctuations. With this type of 
information, assuming that the only change during the occurrence 
of Ca 2+ signals can be captured by a change in the probability that 
a dye molecule be bound to Ca 2+ , numerically simulated images 
of signals with realistic noise can be generated. As mentioned 
before, these images can then be used to determine the minimum 
Ca 2+ current that generates an observable change in fluorescence 
or, as done in this paper, to compute and compare the signal-to- 
noise ratio of different experimental settings. This last type 
application, in turn, can be used to optimize the dye and EGTA 
concentrations with which to observe the signals. Our method is 
particularly useful when experimentalists need to change the dye 
they have been using because it has been discontinued or because 
they want to try a new more efficient one. By means of our method 
experimentalists can compare the detectability properties of 
experiments performed with the old and with the new dye (e.g., 
Fluo-4 and Fluo-8) and also to estimate in which ways the dye 
concentration and/ or illumination intensities could be modified to 
obtain comparable signal-to-noise ratios with both of them. 

In this paper we present a general description of our method 
and then discuss it in more detail when applied to a particular 
example. For the example we obtain microscopy images, under 
stationary conditions (i.e., with no signal evocation), of Xenopus 
lavis's oocytes injected with a single-wavelength Ca 2+ dye (either 
Fluo-4 or Rhod-2) and EGTA. We repeat the experiments varying 
the illumination power or the basal Ca 2+ concentration. 
Analyzing the fluorescence fluctuations in terms of the fluctuation 
model of our method we quantify model parameters. We then use 
the quantified model to estimate the signal-to-noise ratio that can 
be expected from experiments performed using the combinations 
of dye and EGTA that are probed with the stationary experiments. 
We have chosen Fluo-4 and Rhod-2 for the example because they 
have different spectral (one is fluorescein and the other rhoda- 
mine-based) [9, 1 8] and kinetic properties (the on rate of Fluo-4 is 
supposed to be an order of magnitude larger than that of Rhod-2) 
[9,18] and different affinities for Ca 2 + (according to the supplier, 
the dissociation constants are 0.8^Mfor Fluo-4 and 2)iM for 
Rhod-2). Having characterized the detection capabilities of two 
dyes with different spectral properties opens the possibility of using 
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them simultaneously to study different aspects of Ca 2+ signals [19] 
or either one of them in combination with dyes that differentially 
enter the endoplasmic reticulum to monitor cytosolic and luminal 
calcium at the same time [20]. Fluo-4 in combination with EGTA 
has been the dye of choice to observe Ca 2+ puffs in Xenopus Levis 
oocytes [11]. Rhod-2 has been less characterized for this type of 
applications. In particular, the choice of [F/wo4] = 36/iMand 
[EGTA] =90/.(Mhas proven to be adequate to observe Ca 2+ puffs 
[21]. As described later in this paper, using Rhod-2 and EGTA at 
these concentrations, puffs cannot be observed. Our method in 
fact estimates that the expected signal-to-noise ratio of these two 
experimental conditions differs by a factor of two. It places, on the 
other hand, in the same equivalence class in terms of the expected 
signal-to-noise ratio, experiments performed with 
[Fluo4]=36fiMand [EGTA] = 90[iM or [Rhodl] = 90/iMand 
[EGTA] =45{iM, two conditions for which Ca 2 + puffs are readily 
observable. In this example, an analysis of the noisy numerical 
images that can be generated using the quantified fluctuation 
model shed light on the reasons that underlie the differential ability 
of both dyes to detect Ca 2+ puffs at the same concentration. In 
particular, this exploration shows that using the ratio of the 
fluorescence variations during Ca 2+ signals with respect to basal 
fluorescence is not enough to compare experiments performed at 
different dye and/ or EGTA concentrations. We also show how the 
method can be used to quantify the variation of the signal-to-noise 
ratio if the dye and EGTA concentrations, the dye type or the 
illumination power are changed. This particular example then 
highlights the ability of our method as a tool for classification 
purposes and to compare or to improve the detectability 
conditions of different Ca 2 + imaging experimental settings. The 
method, on the other hand, can be used to compare experiments 
performed with different optical set-ups. Other applications of the 
quantified model include the generation of noisy images to 
estimate the smallest detectable signal or the Ca 2 + current that 
underlies a Ca 2 + image. As far as we can tell, this is the first time 
that an analysis of basal fluorescence fluctuations is used for this 
purpose. 

Materials and Methods 

Experiments 

Preparation of Xenopus laevis oocytes. Experiments were 
performed on immature Xenopus lavish oocytes previously treated 
with collagenase and stored in Barth's solution. 

Oocytes were loaded with either Fluo-4 dextran high affinity 
(K d s 800«Af) or Rhod-2 dextran (K d S 2000«M), together with 
an exogenous Ca 2 + buffer, EGTA (Ethylene glycol-bis(2 -ami- 
noethylether)-N,N,N',N'-tetraacetic acid). Intracellular microin- 
jections were performed using a Drummond microinjector. 
Assuming a Ifil cytosolic volume, the final concentration of 
Fluo-4 was 36^A/while the final concentrations of Rhod-2 were 
36/iMor 90/iMand those of EGTA 90/.(Mor 45fiM. We will use 
lower case letters (*, ii, Hi] to distinguish among the three 
combinations of dye type and dye and EGTA concentrations that 
we use in this paper as listed in Table 1. 

We also performed experiments with different (uniform) 
cytosolic Ca 2+ concentrations that were larger than the basal 
value. This was achieved by microinjecting a solution of calcium 
chloride (CaCh). Successive microinjections of this solution were 
applied to a final concentration of 150«Min the oocyte. 

For the observation of Ca 2+ signals, caged IP3 was microinject- 
ed together with the dye and EGTA to a final concentration of 
9(iM. 



Fluo-4, Rhod-2 and caged IP 3 were from Molecular Probes and 
EGTA from Sigma Aldrich. 

Microscopy technique. Experiments were performed using 
an inverted microscope 1X81 connected to a multispectral 
confocal unit Olympus FluoView 1000 in the linescan mode. 
Recordings were made at room temperature. All recordings were 
obtained at the depth of the cortical granules in the animal 
hemisphere of the oocyte focusing with a 60X oil immersion 
objective (NA = 1.35). The dyes Fluo-4 and Rhod-2 were excited 
with the 488«m line of a multiline Argon laser and with the 543«m 
line of a He-Ne laser, respectively. The emitted fluorescence was 
detected in the range of (500 — 60Q)nm for Fluo-4 and 
(555 — 655)nm for Rhod-2 with a PMT detector. Linescan images 
were obtained by scanning along a fixed line (N x = 250pixels) 
within the oocyte (10 — 15)/™. The acquisition rate was fixed at 
lOjU.yper pixel resulting in a scan rate of 3.62m.sper line. Records 
contain a total of N t = 3000lines. The detector was used in the 
photon-counting mode. In this paper we present the results of 4 
types of experiments as summarized in Table 2. 

Observation of IP 3 -mediated Ca 2+ signals. For these 
experiments (which we call Type 0) we use oocytes previously 
microinjected with a Ca 2+ dye, EGTA and caged IP3. To evoke 
the signals, the caged IP3 is photolyzed with a UV pulse (of 
controlled duration and power) using the modification described in 
[21]. This modification allows the entry of ultraviolet illumination 
from a mercury lamp ((350 — 400)nm) while simultaneously 
acquiring fluorescence images with confocal microscopy. 
Throughout the paper we will refer to the "standard illumination 
conditions" as those that are used to excite the fluorescence for the 
observation of IP 3 -mediated Ca 2+ signals. 

Experiments performed under stationary 

conditions. The method that we introduce in this paper is 
based on performing three types of experiments using the dye and 
EGTA concentrations of interest in oocytes where Ca 2+ signals are 
not evoked (i.e., stationary conditions). We call them Type I, Type 
II and Type III experiments (see Table 2). This classification refers 
to the way the experiments are performed. We present as an 
example the application of our methodology to the three 
combinations of dye and dye and EGTA concentrations with 
which we try to observe Ca 2+ signals (sets i, ii, Hi of Table 1). In 
Type I experiments, the fluorescence is collected for 10. 86s along 
a 10.25/.1W1 line using the same illumination power in many regions 
of the same oocyte. In Type II experiments, we select one region in 
one oocyte and obtain several 10.86.? x 10.25,11m linescan images 
each of them for a different illumination power (the "standard" 
one with which Ca 2 + signals are observed and others). To this end, 
the power of the illumination beam is varied over the range that 
the Olympus FV1000 allows with an AOTF. In Type III 
experiments we microinject Cfl 2+ as explained before and obtain 
several 10.86.? x I0.25fim linescan images in different regions of 
the same oocyte and we repeat the data acquisition for different 
final Ca 2+ concentrations. We summarize the characteristics of the 
different types of experiments and list the figures obtained with 
each of them in Table 2. 

Image processing and analyses 

In this paper we present only one figure (Fig. 1) where we 
illustrate the fluorescence distribution obtained in Type 0 
experiments. In this case, the raw fluorescence at each time, tj, 
and point, x,-, along the linescan, F{Xi,tj), is displayed instead of 
the ratio defined in Eq. 1. In Fig. 1, increasing levels of 
fluorescence (increasing calcium levels) are represented by 
increasingly warmer colors. 



PLOS ONE I www.plosone.org 



3 



April 2014 I Volume 9 | Issue 4 | e95860 



Equivalence Classes of Ca + Imaging Experiments 



Table 1. Sets of dye and EGTA concentrations. 





Set 


Dye type and concentration 


[EGTA](fiM) 


M 


[Fluo4\ = 36ftM 


90 


(//) 


[Rhod2]=36/iM 


90 


W 


[Rhod2\=90nM 


45 



doi:1 0.1 371 /journal.pone.0095860.t001 



Most of the results of the paper correspond to Type I-III 
experiments. In those cases we present the results in terms of the 
mean and standard deviation of the fluorescence observed in the 
"bright" fringes {bfj of the analyzed records. In Xenopus Levis 
oocytes typical linescan images show horizontal lines that are 
persistently dark which correspond to the cortical granules 
(Fig. 2A). Bright fringes correspond to the cytosol. The application 
of our method in this case requires a pre-processing of the data so 
as to keep solely the information coming from the cytosol. This 
pre-processing, which is not necessary in other cell types, is 
illustrated in Fig. 2. In order to tell apart pixels corresponding to 
either of these two groups we first compute 



N, 



F(xj) = — F(Xj,tj), and then we calculate the mean, (F} hd , 



and standard deviation, (7e of F(Xf) as: 



CO, 



/),</ = 



1 



b,d ' 



\N X . 



N x _ 



(2) 



(3) 



(Tf- 



\ iebf j = l 



(5) 



where N = Ni,-N t and the sum over i goes over the locations of the 
bright fringes (bfi. 

For each type of experiment (I, II and III) performed with the 
same dye type and dye and EGTA concentrations (sets i, ii, Hi in 
the example discussed in this paper) we typically obtain between 
21 and 88 linescans depending on the experiment type. We then 
put on the same plot the values <r^ vs (F} that correspond to the 
same type of experiment (I, II and III) applied to the same 
combination of dye and EGTA (e.g., i, ii, Hi). Thus, for each set 
{i, ii, Hi) we end up having three plots of ffj, vs (F}. The 
implications and interpretation of this pooling of the data are 
discussed later. We fit the three plots of ff^ vs <_F> with 
polynomials of degree one or two depending on the type of 
experiment using MATLAB's cftool toolbox (The MathWorks, 
Natick, MA). When fitting a curve, this tool gives confidence 
intervals for the various fitting parameters. For the sake of 
simplicity, in the main body of the paper, we only quote the mean 
values obtained. The confidence intervals are listed in the 
Supporting Information (see text SI) file that accompanies the 
paper. 



We finally identify the location of the bright spatial lines, Xi, as 
those for which F(Xj)> (F~) bd — 0.5o> 4i( (rndicated in white in 
Fig. 2B). Once the bright fringes are identified, we only work with 
the pixels of the image whose spatial coordinates correspond to 
those of the bright fringes (Fig. 2C). 

For each linescan image obtained with Type I, Type II and 
Type III experiments we compute the mean, <(-F), and standard 
deviation, c^., of the fluorescence of the pixels in the bright fringes: 



1 N > 



N 



iebfj=l 



(4) 



Fluorescence fluctuations model 

Here we introduce the model with which we describe the 
fluctuations of the fluorescence collected at each pixel of a Ca 2+ 
image obtained using a single-wavelength Ca 2+ dye. With this 
model we analyze the fluorescence distributions obtained with 
Type I, Type II and Type III experiments (Figs. 3-5). By 
determining some key model parameters from fits to these 
experimental data we expect to separate three sources of 
fluorescence fluctuations: (1) variations in the number of dye 
molecules whose fluorescence is collected at each pixel, (2) changes 
in the fraction of such molecules that are Ca 2+ -bound; (3) 
fluctuations in the number of detected photons. We assume that, 
when Ca 2+ signals are evoked (Type 0), the only quantity that 



Table 2. Type of experiments. 







Experiment 


Region of oocyte 


Illumination power 


[c« 2+ U,, 


Caged IP 3 


Figures 


Type 0 




fixed (standard) 


no 


yes 


1 


Type 1 


variable 


fixed (standard) 


no 


no 


3A,4A,5A 


Type II 


fixed 


variable 


no 


no 


3B,4B,5B 


Type III 


variable 


fixed (standard) 


yes 


no 


3C,5C 
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Figure 1. Typical row linescan images obtained in oocytes with EGTA and Fluo-4 or Rhod-2 subjected to the same uncaging 
conditions. (A) For [FluoA] =36fiM, [EGTA]=^iM, (B) for [Rhod2] = 36;iMand [EGTA] =90/iAf, (C) for [Rhodl] = 90/<Afand [EGTA] =45fiM. 
The horizontal and vertical axes correspond to time and space, respectively. The color bar represents the fluorescence intensity (F). The white line 
marks the UV flash. In (A) and (C) several puffs are distinguishable and none can be observed in (B). 
doi:1 0.1 371 /journal.pone.0095860.g001 



changes with respect to the experiments performed under 
stationary conditions is the probability, p, that a dye molecule 
be bound to Ca 2+ at each time and spatial point of the sample. 
We suggest that this probability can either be estimated roughly or 
computed via numerical simulations of the intracellular Ca 2+ 
dynamics in the presence of different Ca 2+ sources. Inserting in the 
quantitated fluctuation model the p(x,t) estimated for a given 
Ca 2+ source the expected signal-to-noise ratio of the correspond- 
ing image can be computed. Repeating the approach for different 
combinations of dye and EGTA (in the example discussed in this 
paper, sets i, it, Hi) allows a direct comparison of the detectability 
properties of different settings in terms of their signal-to-noise 
ratio. 

To build the fluctuation model, we extend the Number and 
Brightness theory [17] taking into account that, for single 
wavelength dyes, both the Ca 2+ -bound and the Ca 2 + -free dye 
molecules emit in the same range of wavelengths (albeit with a 
different intensity) and that the detector introduces an amplifica- 
tion factor. We consider that the fluorescence that is collected at 
each pixel of the linescan image is a random variable that depends 
on the pixel time and on the sources of randomness that we have 
already mentioned. To be more specific, let us call N the random 



variable that represents the total number of dye molecules that 
contribute to the fluorescence at one pixel. The value that N takes 
on at each pixel can be considered as a realization of the random 
variable. As done in [17] we assume that N obeys Poisson 
statistics. Let us call NcaD the number of dye molecules that are 
bound to Ca 2+ at each pixel. This is also a random variable. We 
will call p the probability that a dye molecule is bound to Ca 2+ . In 
general p is space and time-dependent. For the experiments 
performed under stationary conditions we assume that /lis constant 
and spatially uniform. This assumption is dropped for images of 
Ca 2+ signals as explained later. In either case we assume that, 
given N, Nc u d follows a binomial distribution, i.e. 
NcaD = Binomial(p,N). Let us assume that we can describe the 
emitted fluorescence with one molecular brightness for the Ca 2+ - 
bound form of the dye and another one for its Ca 2+ -free form so 
that the total number of photons that reach the detector during the 
acquisition time to eventually give the fluorescence intensity at the 
pixel of interest is given by: No = q\NcaD + l2(N — Ncud)- <7iand 
<72 represent the number of photons per emitting molecule that 
reach the detector during the acquisition time for the Ca 2 + -bound 
and Ca 2+ -free forms of the dye. q\ > q2 and both are increasing 
functions of the power of the laser. 




Figure 2. Image pre-processing applied to Type I —1 1 1 experiments performed in Xenopus leevis oocytes. (A) Typical linescan image 
obtained with a Type I experiment where dark (cortical granules) and bright fringes (cytosol) are distinguishable. (B) Bright fringes (in white) and dark 
ones (in black) of the image in (A) identified as explained in Materials and Methods Section. (C) Final image once the dark fringes have been removed. 
Once this is done we work with all the pixels of the image without distinguishing their time or spatial coordinates. This pre-processing might not be 
necessary in other cell types. The color bar represents the fluorescence intensity {F) both for (A) and (C). 
doi:1 0.1 371 /journal.pone.0095860.g002 
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Figure 3. Fluorescence fluctuations obtained from Type ■ — 1 1 1 experiments performed in oocytes with the set of concentrations (/). The 

mean fluorescence «F» and variance (a 2 F ) are computed as explained in Materials and Methods. The experimental data (black squares) and their 
corresponding fits (black line) are shown for: (A) 84 images obtained in Type I experiments, fit: o-| = 4.8<,F> + 7 '.1; (B) 21 images obtained in Type II 
experiments, fit: cr 2 F = 0.04(F} 2 + 5{F >; (C) 88 images obtained in Type III experiments, fit: a 2 F = 7.27(f) -3 1.4. 
doi:1 0.1 371 /journal.pone.0095860.g003 



Summarizing, we assume that the number of photons that reach 
the detector from the region of the sample associated to each pixel 
can be written as: 



constant of proportionality y, that represents the amplification 
factor, i.e., given JVo, it is: F = yPoisson(No). Combining this 
relation with Eq. 6 we obtain: 



No = (qi — q 2 )Binomial(p,N) + q 2 N (6) 

with TVa Poisson-distributed random variable. We assume that 
<?2/<7l is fixed for each dye and is of the order of the ratio of the 
quantum efficiencies for the Ca 2+ free and the Ca 2+ -bound forms 
of the dye molecule. From Eq. 6 we obtain: 
<N 0 } = {{q 1 ~q 2 )p + q 2 )<Ny and 

r 4/ 0 = ((,(<ll- c l2)P + c l2) 2 +p(l-p)(qi-q2) 2 ^G 2 N for the mean 
and variance of A>o, respectively. Given that N is a Poisson- 
distributed random variable it is a 2 v = (N}. We then conclude 
that a 2 No = (i(q i -q 2 )p + q 2 ) 2 +p(l-p)(qi-q2) 2 )<Ny. Now, 

the detector amplifies the signal and, at the same time, introduces 
some additional noise. In particular, we will assume that, if No 
photons arrive in the detector, the fluorescence intensity that is 
reported at the pixel of interest is a random variable, F, that is 
proportional to a Poisson distributed variable of mean No, with 



F = yPoisson{{q\ — q 2 )Binomial(p,N) + q 2 N) (7) 

with N a Poisson distributed random variable. 

Under these assumptions the mean and the variance (or 

standard deviation squared) of the fluorescence reported at the 
pixel are given by: 

<F} = y((q l -q 2 )p + q 2 )(.N} (8) 

4=r 2 (<^o>+4 0 ) (9) 

= y 2 ((<?! - qi)p +qi+ (<7i - q 2 2 )p + q 2 2 )<N> 




Figure 4. Fluorescence fluctuations obtained from Type I — II experiments performed in oocytes with the set of concentrations (//). 
Similar to Fig. 3 but for set (//). The experimental data and their corresponding fits are shown for: (A) 84 images obtained in Type I experiments, fit: 
<r| = 5.8<F>-4.8; (B) 21 images obtained in Type II experiments, fit: <t| = 0.01<,F> 2 + 6.21<.F >-4.34. In this case the results derived from Type III 
experiments are not shown because no change in fluorescence was observed upon Ca 2+ microinjection for this set of concentrations. 
doi:1 0.1 371 /journal.pone.0095860.g004 
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ABC 




<F> <F> <F> 

Figure 5. Fluorescence fluctuations obtained from Type ■ — 1 1 experiments performed in oocytes with the set of concentrations (///). 
Similar to Fig. 3 but for set (///). The experimental data and their corresponding fits are shown for: (A) 84 images obtained in Type I experiments, fit: 
<r| = 6.2<F> + 1.2; (B) 21 images obtained from Type II experiments, fit: a 2 F = 0.03<f > 2 + 6.02<F>- 15.07; (C) 55 images obtained from Type III 
experiments, fit: 4 = 8.18<.F>-85.31. 
doi:1 0.1 371 /journal.pone.0095860.g005 



Eqs. 8 and 9 can be combined in different ways. Obtaining (N > 
as a function of <-F>, q\ and /from Eq. 8 and replacing it in Eq. 9 
we obtain: 



4 = 7(1 + 



(gl — + 

(q\-qi)p+q2 



(10) 



Using Eq. 8 to write q\ as a function of <.F>, <W> and pand 
inserting it in Eq. 9 we obtain: 



<N) (p + q 2 /q x 



(11) 



Finally, using Eq. 8 to write pas a function of <-F>, <JV> and q\ 
and inserting it in Eq. 9 we obtain: 



a 2 F =y(q\+q2 + i)<F-y-y 2 q iq2 (N-y 



(12) 



We will use either one of these last equations to interpret and fit 
the 4 vs (F)> plots obtained with Type I — III experiments. The 
expression that is used in each case depends on which quantities 
can be assumed to have the same values for the data points that 
are pooled together in each type of experiment (see Table 3). As 
explained in more detail in the Results Section, the combination of 
the fits to the three types of experiments allows the quantification 
of y, q\ at the standard illumination power and<A>. Once we 
have q\, we use the ratio of quantum efficiencies of the Ca 2+ free 
and the Cfl 2+ -bound dye molecules to estimate q 2 /qismd 
determine q 2 . We expect y to depend only on the way in which 
the fluorescence is detected to generate the image, q\ and q 2 to 
depend on the quantum efficiency of the dye and on the laser 
power and (N~) to depend on the dye concentration. Actually, we 
allow for <JV> to vary among different regions of the same oocyte. 
Therefore, with the method we estimate a range of possible <iV> 
values. Now, y, q\ and q 2 at the standard illumination power and 
(Ny should remain approximately the same for any experiment 
(including those of Type 0) that is performed using the same 



experimental set-up, the same dye, in the same concentration and 
with the same illuminating laser power as those performed under 
stationary conditions (Type I III) to estimate them. Thus, we only 
need the fraction of Ca 2+ -bound dye molecules at each spatial 
point and time during a Ca 2+ signal to obtain the expected signal- 
to-noise ratio for each experimental setting. As explained in what 
follows this fraction can be obtained by means of numerical 
simulations of the intracellular Ca 2+ dynamics. 



Numerical simulations 

Determination of the Ca 2+ -bound dye distribution in the 
presence of a localized Ca 1+ source. Numerical simulations 
of the calcium dynamics in the cytosolic medium are performed 
solving a set of reaction-diffusion equations in a spherical volume, 
assuming spherical symmetry, for the following species: Ca 2 + , an 
immobile endogenous buffer (S), a cytosolic Ca 2+ indicator (D) 
and an exogenous mobile buffer (EGTA). In some cases we also 
consider an additional (mobile) buffer (M). A point source of 
calcium located at the origin and pumps (P) that remove calcium 
uniformly in space are also included. The source represents a 
cluster of IP 3 Rs. 

We consider that a single Ca 2+ ion binds to a single buffer or 
dye molecule according to: 



Ca 2+ +X [CaX] 

k off-X 



(13) 



where X represents D, EGTA, S or M, and fc OB _xand k 0 ff-x&re 
the forward and backward binding rate constants of the 
corresponding reaction, respectively. We assume that the total 
concentrations of dye, EGTA, mobile and immobile buffer 
remain constant ([D] T ,[EGTA] T ,[M] T , and [S] T , respectively) and 
that the diffusion coefficient of their free and Ca 2+ bound forms 
are equal. Therefore we calculate the free concentrations, [D], 
[EGTA], [M] and [S] by subtracting the concentration of their 
Ca 2+ bound forms to their total concentrations. Given these 
assumptions, the set of reaction-diffusion equations reads: 
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Table 3. Behavior of the fluctuation model parameters in the experiments performed 


under stationary conditions. 




Experiment 


(N) 


?i 


Pb 


Type 1 


variable 


fixed 




Type II 


fixed 


variable 


fixed 


Type III 


variable 


fixed 


variable 



doi:1 0.1 371 /journal.pone.0095860.t003 



d[Ca 2 +' 



--D C aV 2 [Ca 2+ l 



a" J - ^2 Rcax 

X = D,S,M,EGTA 

\Ca 2+ } 2 

+ source- 5(r) — [P] T - 



lCa 2 +} 2 +k o(r - T 



d[CaD] 
8t 

d[CaEGTA] 
It 



D dye V 2 [CaD}+R c 



= Degta V 2 [CaEGTA] + RquEgta 



d\CaM] 7r 

J --D M V 2 [CaM] + R C aM 



dt 



d[CaS] 



ct 



-Rca, 



(14.a) 

(14.b) 
(14.c) 

(14.d) 

(14.e) 



where D c a , Ddye, ^M an d Degta ar e the diffusion coefficients of 
Ca 2+ ,D, Mand EGTA, respectively. The reaction terms, 
RcaX,^K derived from the kinetic scheme, Eq. 13. 



Rc„x = k on -x[Ca 2+ ]([X] T -[CaX])-k off - X [CaX] (15) 

We assume no flux boundary conditions at r = 20. 5 [im with r the 
radial coordinate. For the source we assume that it consists of nc 
channels that open simultaneously at / = 0 each of which closes 
after a time that is drawn from an exponential distribution with 
mean t open 22]. For the initial condition, we assume that all 
concentrations are homogeneously distributed, Ca 2+ is at basal 
concentration and all species are in equilibrium among themselves 
( J Rc a A' = OforallZ(r = 0)). 

The reaction-diffusion equations are solved using a backward 
Euler method in time and an explicit finite-difference formula in 
space with a 2 nd order expression (first neighbors) for the 
Laplacian. The spatial grid size is dr = 0.041 fim and the time 
step dt=l0fis. The values of the parameters (taken from [18]) 
used in the simulations are listed in Table 4. To compare with 
experimental confocal signals we calculate a weighted average of 
[Cfl-D] along the linescan r = (x,0,0) according to the confocal 
microscope point spread function (PSF). In this case, following 
[23], the confocal signal as a function of x is: 



[CaD](x,t) 

_ J [CaD](r',t) x exp(-2((x-x') 2 +>>' 2 )/uf -2z a /n^dx'dy'dz' ( 16 ) 



where w r = 0.23/im, w z =\A5jxm and V = r w :y^J ■ This is 

the blurred version of the Ca 2+ -bound dye concentration. Having 
[CaD](x,t) we can compute the space and time dependent 
probability that a dye molecule is bound to Ca 2+ during a 



Ca 2+ signal as p- 



[CaD](x,t) 
[D\ T 



Simulated fluorescence distribution 

We generate noisy simulated images using the model described 
in the previous Subsection. Namely, for each position, x, and time, 
t, of the simulation, we draw a stochastic variable, N(x,t), from a 
Poisson distribution whose mean is given by the value, <JV>, 
estimated from the experiments. We then compute the fluores- 
cence at each pixel, F(x,i), following Eq. 7, withA r = N(x,t), 

p= - — - — - - (with [CaD\(x,t) obtained from the simulations as 
[D\ T 

defined in Eq. 16) and the values of y, q\ and <?2 estimated from 
the experiments. 

For each noisy simulated image we compute the signal-to-noise 
ratio as: 



SN-- 



AF 



-signal ' 



-Fba 



(17) 



with F S i sna i and Fbaml the maximum and minimum fluorescence 
values of the simulated image, respectively, and <jj ; computed 



using Eq. 9 with p=pb z - 



[CaD\ b 



An estimate of the signal-to- 



noise ratio can also be obtained replacing F s j gna i and Fb asa l in Eq. 

1 7 by the mean values given by Eq. 8 with p s = - — |' Hj; " a/ and/)/,, 

[D\ T 

respectively. The value of p s can in turn be computed using the 
maximum expected value of the Ca 2+ -bound dye concentration 
during the signal. In this way, the signal-to-noise ratio can be 
written as: 



SN=^/ qi (N) 



Ps-Pb 



(18) 



</2 



neglecting terms of the order of £ = — in front of 1 . Eq. 18 can be 
re-written in terms of quantities that are straightforwardly related 
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Table 4. Parameter values used to solve the simulations introduced in Materials and Methods Section. 





Parameter 


Value 


Units 


Free calcium: 


Dec 


220 


limV 1 




0.1 


uM 


Calcium dye Fluo-4 dextran: 


Ddye 


15 


HmV 1 


kon-D 


240 




koff-D 


180 


s- 1 


ID1 T 


36 


HM 


Calcium dye Rhod-2 dextran: 


Ddye 


15 


rimV 1 


kon-D 


70, 85 




koff-D 


130, 170 


s- 1 


[D] T 


36, 90 


LiM 


Exogenous buffer EGTA: 


Degta 


80 


HmV 1 


kon-EGTA 


5 




koff-EGTA 


0.75 


s- 1 


[EGTAh 


45, 90 


HM 


Endogenous immobile buffer: 




400 


uJVTV 


Koff-i 


800 


s- 1 


Mr 


300 


HM 


Slow endogenous mobile buffer: 


Dm 


27 


umV 1 


Kon-M 


20 


u/vrV 


Koff-M 


8.6 


s- 1 


m T 


250 


HM 


Rapid endogenous mobile buffer: 


Dm 


32 


HmV 1 


K-on-M 


500 




Koff-M 


750 


s- 1 


Mr 


10 


HM 


Pump: 


Koff-p 


0.1 


s- 1 


IPh 


0.9 


HM 


Source: 


«<■ 


6 






20 


ms 



doi:1 0.1 371 /joumal.pone.0095860.t004 



to experimentally accessible parameters. In particular, assuming 
that cji is directly proportional to the intensity of the laser, /, we 

can rewrite q\=q\(I sl ) — , with q\(I s t) the value of q\ at I = I SI , 

1st 

the intensity of the laser at the standard illumination power. 
Replacing pi, as a function of [Ca 



2 + 1 



we can write it as 



Pb = 



[Ca 2 +] 



[Cc 



2+1 



J basal 



In this way, we can compute the SN ratio 



,2+1 



as a function of <iV)>, — and [Ca it> asa i- 

1st 



SN=\/ qi (I sl )-(N) 



Ps-\Ca 1+ } ba J[Cc?+} baMll + K d 



(19) 



^(\+q i ){C a ^\ hclslll /{{Ca^\ hasal + K ll ) + c(\ + Ccn) 



Results 



In the Materials and Methods Section we introduced a method 
to derive a quantitative description of the fluorescence fluctuations 
of images obtained using single-wavelength Ca 2 + dyes. In this 
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Section we derive this quantitative description for three situations 
(sets (i)-(iii) of Table 1). The first one corresponds to a choice of 
dye and concentrations that has proven to be adequate for the 
observation of Ca 2+ puffs [21], (set (i)). The other two involve 
using a dye that has been less characterized for the observation of 
this type of signals, Rhod-2 (sets (it) and (tit)). We have chosen 
these three sets to show that the fluctuation model can be used to 
classify Ca 2+ imaging experimental settings and that the classifi- 
cation is able to discriminate between experimental settings with 
clearly different signal-to-noise ratios (when probed with Type 0 
experiments). The differences in the signal-to-noise ratios for these 
three situations are apparent in Fig. 1. We show in this figure 
linescan images obtained in oocytes after photoreleasing IP 3 with 
the same uncaging pulse (the occurrence of which is marked with a 
white line in the images). Fig. 1A corresponds to set (t), Fig. IB to 
set (it) and Fig. 1C to set (tit). Puffs can be observed in Figs. 1A 
and 1C, but not in Fig. IB. In fact, we repeated the experiments 
for the conditions of Fig. IB in many different oocytes and could 
never observe isolated signals using Rhod-2 and EGTA at these 
concentrations. As we show in this Section, our method puts the 
experiments of Figs. 1A and IB in two separate classes. It also 
predicts that the signal-to-noise ratio of the experiment in Fig. 1 C 
should be similar to that of Fig. 1A. More interestingly, the 
numerical simulations show that this separation in two different 
classes is not achieved if one compares the Ca 2 + -bound 
distribution without including the various sources of fluctuations 
that our method discriminates. In this Section we also discuss the 
rationale for analyzing the data of the various experiments in the 
way that we do it. 

Outline of the method to quantitate the contribution of 
the different sources of fluorescence fluctuations in 
images obtained using single-wavelength Ca 2+ dyes 

Before advancing with the particular application of our method 
to sets (i)-(iii), we first oudine how to proceed in a generic 
situation. Let us suppose we want to obtain images of Ca 2+ signals 
using a single-wavelength dye (e.g., we want to perform an 
experiment of Type 0) and that we wish to evaluate a priori the 
performance of the experimental set-up for a certain combination 
of dye, illumination power and dye and EGTA concentrations 
(i.e., for an experimental setting). Then, we microinject the desired 
quantities of dye and EGTA in different oocytes and perform 
Type I, Type II and Type III experiments as described in the 
Materials and Methods Section. For each experimental type we 
process the data as described before and obtain plots of 
fluorescence variance, (T F , as a function of mean fluorescence, 
(Fy. For classification purposes a simple possibility is to compare 
the three curves obtained for the dye of interest with those 
obtained for another dye or at other concentrations at which it is 
known that Ca 2+ puffs can be observed. Experiments performed 
with different set-ups can also be compared. For a more 
quantitative comparison we use one among Eqs. 10-12 to fit 
each of the three ff^vs<-F> curves obtained experimentally. We 
discuss in the next Subsection what equations are applicable in 
each case. The three fits should serve to quantify the parameters y, 
q\ and qi (at the standard illumination power) and to determine a 
range of values for the mean number of dye molecules, (N) for 
the experimental settings that are being probed. When using the 
same experimental set-up and the same combination of dye type, 
dye and EGTA concentrations, these four parameters should 
remain approximately the same for Type 0 and for Type I-III 
experiments. The signal-to-noise ratio during a Ca 2+ signal (i.e., in 



a Type 0 experiment) depends on these four parameters and on 

Ps= - — |'"^" g/ ( see Eq. 18). Thus, once the four parameters are 
[D\ T 

determined with Type I-III experiments, it is only necessary to 
have an estimate of the fraction of dye molecules that are bound to 
Ca 2+ at the location of the signal maximum in order to compute 
the signal-to-noise ratio that can be expected in a Type 0 
experiment. This fraction can be estimated roughly or by means of 
numerical simulations, as explained in Materials and Methods 
Section. This outline shows that the application of our method 
allows an a priori (quantitative) estimate of the signal-to-noise-ratio 
for a given experimental setting. We now describe in detail the 
application of this approach to the sets (i)-(iii) explored in Fig. 1 . 

Fluctuation analysis of experiments with 

[FluoA] = 36^Mand [EGTA] = 90/iM 

We show in Figs. 3 A, B and C the variance of the fluorescence, 
op, as a function of its mean, derived from Type I-III 

experiments performed on oocytes microinjected with the set of 
concentrations (t) of Table 1 . 

Fig. 3A collects the data obtained in several regions of the same 
oocyte with the standard illumination power and without 
Ca 2+ microinjection (Type I see Tables 2 and 3). Thus, q\ and 
q2 take on the same values for all the data points of Fig. 3A. The 
fraction of Ca 2+ -bound dye molecules, p, could vary from region 
to region due to a non-uniform distribution of basal cytosolic 
Ca 2+ . We observe in Fig. 3A, however, that a 2 F and <-F> are 
linearly related and that the op vs(Fy curve practically goes 
through the origin (the ordinate to the origin given by the fit is an 
order of magnitude smaller than the values of a 2 ? obtained in the 
experiments). Given Eq. 10, this is expected if variations in the 
basal value of p (pi) in different regions of the same oocyte do not 
have a noticeable effect on the ratio op/(Fy. We then use Eq. 10 
to describe the curve displayed in Fig. 3A and equate 

. Y j + (gi -q 2 )Pb+q 2 \ to the sl ofthe linear fit) u 4 g This 
V (q\-q2)pb+q2j 

is the first piece of quantitative information that we derive from the 
experiments. It is implicit in this interpretation of the data that 
different points along the <sp vs <-F)curve have different values of 
<JV> which is consistent with having slightly different dye 
concentrations in different cytosolic regions of the same oocyte 
(Table 3). 

We show in Fig. 3B the plot of a 2 F vs <-F> derived from the 
analysis of linescan images obtained in Type II experiments. In 
this case, since we do not pool together data coming from different 
regions, we assume that p=Ph and (N) take on approximately 
the same value for all data points (see Table 3). Regarding q\and 
q2 we can only assume that their ratio is constant but not their 
individual values since we are varying the laser power. We then 
use Eq. 1 1 to interpret the fits to the data. In particular, fitting the 
data points with a second degree polynomial we derive y = 5. 
Comparing these estimates with the slope derived from the fit of 

Fig. 3 A we conclude that y( 1 + { £lZ3^E!LL3l \ ~ y for the 

experiments of Fig. 3A. Thus Eq. 10 reduces to: 

4 (2°) 

under basal conditions for the standard illumination power used in 
Fig. 3A. The estimates of the amplification factor y = 5 also agree 
with an analysis of the detectors behavior in the photon counting 
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mode and in the absence of a fluorescence stimulus (data not 
shown). In particular, after subtracting the zero counts, the mean 
number of counts reported by the detectors in darkness is 
consistent with this amplification factor, namely, it is 
(,Fdarky\pj:0^-6- If the zero counts are included as well we find 
<^<&mt)0.61 which shows that the detectors do not introduce a 
noticeable offset in the reported fluorescence. 

We show in Fig. 3C the plot of cr^vs (Fy derived from the 
analysis of linescan images obtained with experiments performed 
using the same (standard) illumination power as in Fig. 3A but in 
oocytes microinjected with dye, EGTA, and different amounts of 
Ca 2 + to change the cytosolic Ca 2+ concentration, [Ca 2 + ] cyt . 
Data coming from the same batch of oocytes but with no added 
Ca 2 + are also included in the plot. In Fig. 3C not all the data 

[CaD] 

points correspond to the same value oip = -j— — . Furthermore, 

[D\ T 

the value of <JV> cannot be assumed to be the same for all the 
points either (see Table 3). In any case, if (Fy is large enough so 
that y(l+qi+q 2 )(F)> >y 2 q\qi(Ny, Eq. 12 reduces to: 

4 = V(.11+12 + 1)<F> (21) 

A similar relationship is obtained from Eq. 10 if p is large 
enough so that (q\ — qi)p> > qi, or equivalently, p > >qi/q\ (for 
situations like those analyzed in this paper in which q\> > qi). 
Thus, using Eq. 21 to fit the points with large enough <F> of 
Fig. 3C we estimate y(l +<7l +<72)- We obtain 
y(l +qi + qi) «7.27 — 7.54 depending on how many points we 
keep to do the fit «F> > 12 as illustrated in Fig. 3C or <F> > 18, 
respectively). Using the values of y estimated from Fig. 3B and 
assuming that qi/qi < < 1 we derive q\ from these fits. We obtain 
qi =0.45 — 0.5. Assuming that qi/q\ is of the order of the ratio of 
quantum efficiencies estimated for Fluo-3 in [24], qijq\ =0.025, 
we can derive qi as well. We obtain qi = 0.01 1 for;?] = 0.45. These 
values correspond to the standard illumination power since they 
are derived from Type III experiments (see Table 2). 

From Fig. 3A a range of possible mean values, <A>, can be 
estimated using Eq. 8 with the previously inferred values, y = 5, 
qi =0.45 and qi/q\ =0.025. In particular, using pi, = 0.125 (which 
corresponds to [Cc? + ] basaI = 0.1(iM and Kj^O.SfiM) we find 
that the range of mean fluorescence values of Fig. 3A 
(14<<i 7 ><22) corresponds to 41<<iV><65. If we consider 
Pb = 0.03, instead, (which corresponds to [Ca 2+ ] haml = 0.02/iM) 
we obtain 1 13 < <JV> < 177. 

Summarizing, from the fits analyzed so far we have estimated 
y = 5, q\ =0.45, i/2 = 0.01 1 and a range of <JV> for different values 
of pb for set (i). 

Rough estimate of the expected signal-to-noise ratio for 
the observation of Ca 2+ puffs for experiments with 
[FluoA] = 36/^Mand [EGTA] = 90/iM 

Having obtained y, q\, qi and (N) we are in a position of 
estimating the signal-to-noise ratio as a function of the fraction of 
Ca 2 + -bound dye, p, during a signal. Here we only give a rough 
estimate. A more accurate description can be obtained with 
numerical simulations as described later. To this end we use Eq. 
18. Given the previous estimates of the model parameters we 
obtain AF /c basal W{Ps-Pb) for p A = 0.125 and <iV> = 45. Thus, 

already at (p s —pb)6.2 it is Irsbasa?' ant ^ we ex P ect tne sig na l to 
start to be distinguishable. An estimate of the underlying 



Ca + source that produces a difference />., — p/,0.2 can be obtained 
with numerical simulations. 

Fluctuation analysis of experiments performed with 
Rhod-2 

We now repeat the experiments and analyses described in the 
previous Subsections but for oocytes microinjected with two sets of 
Rhod-2 and EGTA concentrations: set (it) with the same dye and 
EGTA concentrations as in set (i) and with set (Hi) with a larger 
dye and smaller EGTA concentrations. Given that the same laser 
and the same standard illumination conditions are used for sets (it) 
and (Hi), we expect that the fluctuation model should be 
characterized by the same values of y, q\ and qi at the standard 
illumination power. The value of (N~), however, should be 
different due to the different values of [Rhod2] . 

We show plots of the variance a 2 F , as a function of the mean 
fluorescence, <-F>, derived from Type I II experiments performed 
using set (it) in Figs. 4A and 4B and from Type I — III experiments 
performed using set (Hi) in Figs. 5A, 5B and 5C. 

We show in Figs. 4A and 5A the plots derived from Type I 
experiments. We observe that for both sets of concentrations <s 2 F 
and (Fy are linearly related and that the a^vs^Fy curve goes 
through the origin. Furthermore, for both types of experiments the 
fit to the data approximately gives the same slope (~6). Thus, 

y ( 1 + —J — _2_Pj' I2 j ~ g r or typgg of experiments which is 

V (?i-«2)p»+9iy 

a consequence, as in the case of Fluo-4, of the fact that the second 
term inside the parenthesis is much smaller than one. 

We show in Figs. 4B and 5B data points obtained from Type II 
experiments. As done for Fluo-4 we assume that p=Pb and (Ny 
take on approximately the same value for all the points of each 
figure and that the ratio, qi/q\, but not the individual values, q\ 
and qi, is constant. In this case we observe a fundamental 
difference between the plots of Figs 4B and 5B. Namely, while the 
G F vs(Fy curve has a parabolic shape in Fig. 5B, it remains a 
straight line in Fig. 4B. This latter behavior is different to the one 
observed for Fluo-4 too (Fig. 3B). Using either a linear or second 
degree polynomial to fit the points of Fig. 4B and interpreting the 
results with Eq. 1 1 we derive y = 6. We obtain a similar value from 
the fit (in this case, a second degree polynomial) of the points in 
Fig. 5B. We then conclude that for the experiments performed 
with Rhod-2 (that use the He-Ne laser) it is y = 6. The different 
behavior observed in Figs 4B and 5B can then be attributed to the 
fact that the nonlinear part of the a 2 F vs <,F> curve is not reached 
for the laser powers probed in case (it). We discuss later possible 
causes for this different behavior which can be used to discard 
combinations of dye and EGTA concentrations that are not good 
for the visualization of localized Ca 2 + signals. 

We show in Fig. 5C the data points obtained from experiments 
of Type III. We do not show the corresponding figure for case (it) 
because we did not observe any variation of the fluorescence upon 
Ca 2+ injection. The reason that underlies the different behaviors 
observed in Figs. 4B and 5B may also underlie the lack of 
fluorescence variations observed upon Ca 2+ microinjection for 
case (it). From the slope of a linear fit to the points with large 
enough <.F> of Fig. 5C we estimate y(q\ +qi+ 1) as before. Using 
y = 6 and assuming that qi/q\ < < 1 we obtain q\ =0.36 at the 
standard illumination power. 

We now estimate the range of possible mean values, (N), from 
the data of Figs. 4A and 5A. To this end, we use Eq. 8 with y = 6, 
^1=0.36 and qijq\= 0.07. In particular, for case (it) from Fig. 4A, 
we obtain 39<<W><98 using p=Pb = 0.05 (which corresponds 
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to [Ca 2+ ] hasal = 0ApM) and 58<<JV><145 using #, = 0.01 
(which corresponds to [Ca 2 + ] baml = 0.02/.iM). These values are 
consistent with the ones derived for Fluo-4 at the same 
concentration of dye and EGTA (Fig. 3A). For case (tit) we 
obtain, from Fig. 5A, 98<<7V><177 for /7 fe = 0.05and 
145<<JV><260 forp/, = 0.01. These values of <#> are between 
1.8 and 2.5 times larger in Fig. 5 A than in Fig. 4A. 

We now estimate the signal-to-noise ratio as a function of the 
fraction of Ca 2+ -bound dye,/).,, during a signal (i.e., the expected 
value during a Type 0 experiment). We again use Eq. 18. Setting 
<iV>=115and pi,~0.05 for set (lit) we obtain /oho* i~^(P~ Pi>)- 
Given the different kinetic properties and dissociation constant 
between Fluo-4 and Rliod-2 we do not know a priori by how much 
the Ca 2+ -bound dye concentration in the presence of the same 
source would differ for set (i) and for sets (ii) or (tit). This can be 
estimated by means of numerical simulations. In particular, 
numerical simulations of the intracellular Ca 2+ dynamics when 
there is a localized Ca 2+ source performed as explained in 
Materials and Methods Section show that the ratio 
p(Rhod2)/p(Fluo4) is slightly larger at the peak of the signal 
than at basal conditions when simulations for sets (i) and (Hi) are 
compared (see next Subsection). In particular, for the simulations 
displayed in Figs. 6A and 6C it is p s (Rhod2)/p s (Fluo4)0.5 while 
Pb(Rhodl) / ph{FluoA) = 0.5/0.125 = 0.4. Using these estimates we 
conclude that the source that would give a signal-to-noise ratio of 
the order of 2 for Type 0 experiments performed with set (i) 
(p s (Fluo4)~0.32) corresponds to p s (Rhod2)~0A6 (for Type 0,set 
(Hi)) experiments. Inserting />, = 0.16 and />/,~0.05 into 
^loh ai ~^7(Ps ~ Pb) we determine that such a source would give 
a signal-to-noise ratio of the order of 1 .9 for Type 0 experiments 
performed with set (Hi). 

This estimate drops by a factor of 1.6 if (,N~> is decreased to the 
values that can be expected for the set (ii). Thus, based on the 
rough estimates of the signal-to-noise ratios derived so far, the 
application of our method puts the experiments of case (i) in the 
same equivalence class as those of case (Hi). So, in principle, 
localized Ca 2+ signals should be equally detectable with either one 
of these experimental settings. The smallest concentration of 
EGTA used in the latter could be compensating the larger 
dissociation constant of Rhod-2 and be the cause for this to 
happen. In fact, we prove that this is the case with numerical 
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simulations that allow us to estimate p for a very small signal. 
According to these rough estimates the set of concentrations (ii) 
belongs to a different class. 

Numerical simulations of the model to compute signal- 
to-noise ratios of Type 0 experiments and establish 
equivalence classes 

We now show the results of using numerical simulations to 
compute the expected signal-to-noise ratio for the experimental 
conditions explored in Fig. 1. For this particular example we 
simulated the Ca 2+ , dye and buffer dynamics as described in 
Materials and Methods Section using the parameters of Table 4 
but witiiout including the mobile endogenous buffers listed in that 
Table. We used the same Ca 2+ source in the three cases. In 
particular, we chose a situation with n c = 6 IP^Rs that were 
initially open and subsequently closed randomly with mean time, 
t 0 pen = 20ms. We have chosen this source because it gives signal 
amplitudes of the order of the smallest detectable one (SN2). We 
explored other sources obtaining similar results as those illustrated 
in this Section (data not shown). We show in Fig. 6 the blurred 
Ca 2+ -bound dye distributions obtained with the simulations 
where (A) corresponds to set (i), (B) to set (ii) and (C) to set (Hi). 
Once we had these noiseless plots, we replaced the value of the 
Ca 2+ -bound dye concentration at each "pixel" by a random 
fluorescence value as explained in in Materials and Methods 
Section. We show the resulting noisy images in Fig. 7. To obtain 
more realistic images we also added dark fringes that correspond 
to the granules in the oocyte. Figs. 7A and 7D correspond to 
simulations with set (i), Figs. 7B and 7E to set (ii) and Figs. 7C and 
7F to set (Hi). In order to go from the Ca 2+ -bound dye to the 
fluorescence distribution we used y = 5 and (71 =0.45 for Fluo-4 
and y = 6 and q\ =0.36 for Rhod-2, as derived from the previous 
analyses. For set (Hi) (which has [Dye] = 90pAF), we chose <W> 
among the values estimated from Fig. 4C, namely, (.N) = 100. For 
comparison purposes, we chose <(A f ) = 40»36/90 x 100 for sets 
(i) and (ii) (which have [Dye] = 36/iM). Figs. 7A-C were obtained 
using values of q 2 such that their ratio with respect to (/jwas equal 
to the ratio of quantum efficiencies estimated in [24], i.e., q 2 =0.01 
(q 2 /q l =0.025) for Fluo-4 (set (i)) and ^ = 0.02 (q 2 /qi =0.07) for 
Rhod-2 (sets (ii) and (Hi)). Figs. 7D-F were done using q 2 = 0 to 
analyze the role of the Ca 2+ -free dye fluorescence on the signal- 
to-noise ratio for this particular example (see Discussion). 




0 0.05 0.1 0.15 0 0.05 0.1 0.15 0 0.05 0.1 0.15 

Time [s] Time [s] Time [s] 

Figure 6. Simulated ( 1/ -bound dye concentration during a Type 0 experiment in a linescan image. Blurred Or -bound dye 
concentration obtained for: (A)[/Vuo4] =36/iM, [EGTA] =90/iM(set (/)); (B) [Rhodl] = 36^M, [EGTA] =90^M(set (//)); (C) [Rhodl] = 90/iM, 
[EGTA] =45/iM(set (///)). All other parameters are as given in Table 4. In all the simulations a puff involving the simultaneous opening of 6 IP 3 Rs, 
n c = 6, occurs at time t = 0. 
doi:1 0.1 371 /journal.pone.0095860.g006 
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Comparing Figs. 7A-C we observe that the puff is less 
distinguishable when using Rhod-2 and EGTA at the same 
concentrations as in the simulations with Fluo-4 (Figs. 7A and 7B, 
respectively). The puff becomes distinguishable when the concen- 
tration of Rlrod-2 is increased and that of EGTA is decreased 
(Fig. 7C). These results agree with the experiments of Fig. 1. These 
qualitative observations are quantified in the second column of 
Table 5 where we show the signal-to-noise ratio computed as 
explained in Materials and Methods Section. There we see that 
the signal-to-noise ratio is smallest for the set of concentrations (it) 
and that the one obtained for the set (Hi) is only slightly smaller 
than the one obtained for the set (i) so that a similar level of 
detectability of Ca 2+ signals can be expected for these two 
experimental conditions (in Type 0 experiments). Thus, the 
numerical simulations of the model confirmed our previous rough 
comparison between the experimental settings explored in Fig. 1 
putting in the same equivalence class the conditions of Fig. 1A and 
Fig. 1C and in a different one those of Fig. IB. 

Using the model to choose adequate experimental 
parameters 

Once the model is quantified for a certain dye, standard 
illumination power and dye and EGTA concentrations, it can then 
be used to determine how changes in these experimental 
parameters affect the expected signal-to-noise ratio. As expected, 
it follows from Eq. 18 that the signal-to-noise ratio increases with 
<JV> (i.e., with the dye concentration) and with q\ (i.e., with the 



illumination power or, equivalently, the laser intensity, /) while it 
decreases with pi, (i.e., the concentration of basal Ca 2+ , 



[G 



.2 + 1 



The quantified fluctuation model, however, gives 



more information than that. Namely, it can be used to select the 
experimental conditions in a more quantitative way. As an 
example, we show in Fig. 8 plots of the signal-to-noise ratio, given 
by Eq. 19, as a function of <A> (A), ///„ (B) and [Cfl 2 + ] w (C) 
for two dyes. The aim of this figure is two-fold. First, we illustrate 
how to use our method to select good experimental conditions for 
a given dye (in the example, Fluo-4). The second goal is to 
illustrate how to select good experimental conditions when 
planning the replacement of a dye by another that is excited with 
the same laser but that has different photo-physical properties (in 
the example, Fluo-4 and Fluo-8, respectively). This approach can 
be applied to ensure comparability between experiments per- 
formed with the "old" and the "new" dye. For these plots we first 
computed the signal-to-noise ratio using Eq. 19 with the values 
y = 5, q\(I st ) = 0.45 and ^C/s;) = 0.01 1 , determined from the 
application of our method to set (t) and <iV> = 40 in (B) and (C), 
the standard illumination intensity, I = I st , in (A) and (C) and 



[G 



,2 + 1 



-0.lf.iM in (A) and (B). For these three curves (shown 



with solid lines in the figure) we computed the signal-to-noise ratio 
assuming that the fraction of dye molecules that are bound to 
Ca 2+ at the location of the signal maximum is p s = 0.3 (the 
smallest detectable signal amplitude according to the previous 
discussion in Result Section). The vertical dotted lines indicate the 
range of <A r ) values determined from the application of our 
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Figure 7. Simulated fluorescence during a Type 0 experiment in a linescan image. The fluorescence is computed from the Or 1 -bound 
dye distributions of Fig. 6 as explained in Materials and Methods Section. (A) and (D) correspond to set (i), (B) and (E) to set (ii) and (C) and (F) to set 
(iiii). (A), (B) and (C) are obtained using qi/q\ equal to the ratio of quantum efficiencies estimated in [22]. (D), (E) and (F) are obtained setting qi =0, 
i.e., they do not take the Ca 2+ -free dye fluorescence into account. In all cases the values of q\ are the ones derived from the fluctuation analyses for 
the standard illumination power. 
doi:1 0.1 371 /journal.pone.0095860.g007 
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Table 5. SN ratios obtained in the noisy simulations. 



Set q 2 = 0 92#0 

M 2.7 2.5 

(/;) 2.3 1.5 

(Hi) 3.1 2.1 



doi:1 0.1 371 /journal.pone.0095860.t005 

method to set (i) in Fig. 8A, the standard illumination condition, 

— =1, in Fig. 8B, and, in Fig. 8C, the value of [Ca l+ ] basal we 

have been using for set (t), [Ca 2+ ] ham/ = 0.1 fiM (which corre- 
sponds to pb = 0.1, Kd = Q.%nM). From the observation of the solid 
line curves in Fig. 8 we conclude that moderate increments in I 
would not result in noticeable changes in the signal-to-noise ratio 
of Type 0 experiments performed with Fluo-4. Something similar 
happens with [Ca 2 + ] haml , a parameter whose change would 
directly affect the type of signals that can be evoked. The natural 
variations in (N") that can be encountered by changing the region 
of the oocyte where the experiments are performed, on the other 
hand, can lead to a ~20% improvement in the signal-to-noise- 
ratio. 

The dashed curves in Fig. 8 are the equivalent of the solid line 
ones but for Fluo-8. In order to compute them we used 
Kj = 0A/j.M and assumed that Fluo-8, when bound to Ca 2+ 
and excited with a 488«m wavelength, is approximately twice as 
bright as Fluo-4 while the brightness of the Ca 2+ -free form of both 
dyes is approximately the same (see e.g. Fluo-8 data specification 
sheet, Teflabs). The last two assumptions imply that we replaced 
qi by 2q l (I st )I/I st = Q.9U/I s , and £ = 0.00125 when going from 
Eq. 18 to Eq. 19. We estimated the value of/)., for Fluo-8 simply by 
assuming that it would differ from the value used for Fluo-4 
(p s = 0.3) as if Ca 2+ and the dye were at equilibrium at the peak of 
the signal. Namely, we determined [Ca 2+ \ sigml at the peak from 

Ps (Huo4) = 03 = [Ca 2+ }^ ml /([Ca 2 % tgml +K d (Fluo4)) and 

then computed p s (Fluo8) using the same formula but with 
Kj(FluoS). We obtained p s (Fluo8) sO.5. As in the case of Fluo- 
4 we used (N) = 40 in (B) and (C), the standard illumination 



intensity, 1 = I st , in (A) and (C) and [Ca 2+ ] hasal = 0.lfiM in (A) and 
(B). From a comparison of the solid and dashed curves we can then 
estimate by how much experimentally accessible parameter values 
should be changed in order to obtain equivalent signal-to-noise 
ratios with the old (Fluo-4) and the new dye (Fluo-8). It can also be 
used to determine the expected difference in signal-to-noise ratios 
when using both dyes in the same concentration. 

Discussion and Conclusions 

Optical techniques and Ca 2+ fluorescent dyes offer the 
possibility of observing intracellular Ca 2+ signals with minimum 
disruption. The fluorescence changes that occur during signals are 
related to the Ca 2 + -bound dye concentration, a quantity that not 
only depends on free Ca 2+ but also on the ability of the dye to 
overcome other Ca 2+ trapping mechanisms such as buffering. 
Furthermore, when trying to observe IP 3 -evoked localized 
Ca 2+ signals, exogenous buffers such as EGTA are typically 
introduced in the cell to disrupt the Ca 2+ -mediated communica- 
tion between clusters of channels. This additional buffer also 
competes with the dye for Ca 2+ and can degrade the signal. There 
is great interest in studying local IP 3 R-mediated Ca 2+ signals 
(puffs) since they constitute the building blocks of more global 
signals. Puffs are highly stochastic. In order to have an accurate 
statistical description of their properties it is necessary to have 
experiments that report the occurrence of puffs with similar levels 
of accuracy over the whole range of event sizes. Being able to 
estimate the signal-to-noise ratio that can be expected from an 
experimental setting can at least provide information on the 
reliability of the event size distribution in the region of small 
events. In this paper we have introduced a method that provides 




Figure 8. Dependence of the signal-to-noise ratio on experimentally accessible parameters. Plots of the signal-to-noise ratio using Eq. 1 9 
with the parameter values determined for Fluo-4 (solid line) and for Fluo-8 (dashed line) as a function of the mean number of dye molecules, <iV> (A), 
the normalized laser intensity, I/I sl (B) and [Cct + ] lmal (C). The vertical dotted lines indicate the range of <W> values and the values of I/I st and 
[Ca 1+ } basa i that correspond to set (/) at the standard illumination power. See text for more details. 
doi:1 0.1 371 /journal. pone.0095860.g008 
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an a priori estimate of the signal-to-noise ratio of experiments that 
use single wavelength Ca 2+ dyes and EGTA. These are the dyes 
that are used to observe IP :i R-mediated Ca 2+ signals that are 
evoked via the photorelease of caged IP 3 with UV illumination. 

Our approach to probe the performance of Ca 2 + -imaging 
experiments is based on the Number and Brightness (N&B) 
method of [17,25]. It entails performing a series of experiments 
under stationary conditions (Type I — III) from which the param- 
eters of a model of fluorescence fluctuations can be quantified. 
Combining the quantified fluctuation model with simulations of 
the dynamics of intracellular Ca 2+ in the presence of a localized 
Ca 2+ source, the expected signal-to-noise ratio in signal evoking 
(Type 0) experiments can be estimated and images with realistic 
noise can be generated. In our problem there are additional 
sources of fluctuations with respect to the traditional N&B method, 
even under stationary conditions. Besides variations in the total 
number of fluorescent molecules and in the number of detected 
photons, variations in the fraction of Ca 2 + -bound dye molecules 
also contribute to fluorescence fluctuations at each pixel. 
Furthermore, this fraction depends on the number of free Ca 2+ 
ions in the region of interest which is also a dynamic variable. Our 
method involves a simplification in this sense. Namely, we assume 
that the probability, p, that a dye molecule be bound to Ca 2+ is 
fixed under basal (stationary) conditions and that its dynamics is 
governed by deterministic equations (Eqs. 14) during Ca 2 + signals. 
Thus, in the absence of signals, the fluorescence fluctuation model 
of our method has five unknown parameters, the amplification 
factor of the detector, y, the number of emitted photons per free 
and Ca 2+ -bound dye molecule that reach the detector, (^and qi, 
the probability that a dye molecule be bound to Ca 2+ , p, and the 
mean number of dye molecules that contribute to the fluorescence 
at a pixel, <iV>. These unknown parameters reduce to four if 
qi/li is known a priori, as assumed in this paper. When 
quantifying the model of our method we also assume that p and 
<(iV) can vary from region to region of the sample (in our case, the 
oocyte) and between samples (i.e., between oocytes). This increases 
the number of quantifiable parameters with respect to the N&B 
method. That is why our method involves the performance of 
more experiments than those of N&B. In one of these additional 
experiments the cytosolic [Ca 2+ ] is varied which, in turn, varies/;. 
The purpose of our method is more restricted than that of the 
N&B technique and this simplifies the analyses. We are not trying 
to build a pixel-by-pixel map of [Ca 2+ ]. We are just interested in 
having a realistic quantified model of the distribution of 
fluorescence fluctuations. 

Having a realistic fluorescence fluctuation model can be helpful 
to estimate the Ca 2+ current that underlies an image. One 
possibility is by means of what is called a "forward approach" 
[26] in which the experimentally obtained image is direcdy 
compared with one that is generated via numerical simulations as 
those of Figs. 7A-F. Our fluorescence fluctuation model can also 
be useful for backward methods in which the current is inferred 
directly from the image [27,28]. Even in this type of approach, 
once the current is inferred a numerically simulated image is 
generated for direct comparison with the experiments and this 
requires the addition of noise which is usually done in an ad hoc 
way [29,30]. 

Our method can be used to identify small fluctuations that are 
due to Ca 2+ release through one or a few open channels. This is 
relevant for the use of optical techniques to infer single channel 
kinetics (the so called "optical patch clamping") [31]. Being able to 
identify single channel openings is particularly important in the 



case of IP3RS. According to a variety of experiments IP^Rs diffuse 
on the membrane of the endoplasmic reticulum [32,33] and while 
some observations seem to indicate that IP3R clustering occurs as 
a consequence of stimulation with IP :i [15,16], other Ca 2+ 
imaging experiments indicate that Ca 2+ puffs occur at fixed 
locations in the cell [34]. Being able to tell apart spontaneous 
fluctuations from small changes in fluorescence due to Ca + 
release would certainly be of help to solve this apparent paradox. 
Small changes in fluorescence with respect to the background also 
need to be detected when trying to analyze the dynamics of 
luminal Ca 2+ during the occurrence of localized signals as done in 
[35]. In [35] fluctuations due to Ca 2+ release are distinguished 
from spontaneous fluorescence fluctuations by means of an 
analysis that assumes that fluctuations follow a Poisson distribu- 
tion. A model like the one constructed with our method, that 
quantifies the fluctuations of each specific experimental setting and 
that separates the background noise due to fluctuations in the 
number of fluorescent molecules from those in the number of 
detected photons that is present in all images [36] would certainly 
be a useful tool to identify whether a small fluctuation in luminal 
Ca 2+ is due to Ca 2+ release into the cytosol or not. 

Self-consistency tests of the method 

In this paper we have used our method to compare the signal 
detectability properties of experiments performed with two 
different dyes and at different concentrations. Before discussing 
the signal-to-noise ratios that we could estimate for the sets 
we do first some self-consistency tests of the method. The first test 
consists of determining a range of [Ca 2+ ] haml values that are 
compatible with the fluctuation model parameters that we 
estimated with the method for sets In particular, using 

the values of y and of q\ and 92 a t the standard illumination power 
we can establish a range of possible values of pb (the probability 
that a dye molecule be bound to Ca 2+ under basal conditions) 
that are compatible with the results derived from Type III 
experiments. We do it first for set (i). To this end we use Eq. 12 
with y(l +q\ + q 2 ) «7.3 — 7.5 on the points of Fig. 3C with low 
enough fluorescence to guarantee that they correspond to basal 
Ca 2+ conditions. In particular, we do it for the points with 
(Fy< 12 and with<(.F> < 18. Applying Eq. 12 to these points we 
obtain for each of them a value of the ordinate, — yq\q 2 (Ny. 
Assuming that p=Pb for these points we use Eq. 8 to write 
y 2 qiq2(N} = y(Fyq 2 /(pb + q2/qi)- From this equality we deter- 
mine a set of possible values of p/,. In this case we find: 
p/,x(Q.l — 0.1)q 2 /q l . This means that, for set (t), the contribution 
of the fluorescence of the Ca 2+ -bound dye can be comparable to 
that of the free dye under basal conditions. Using pb = 0. lq 2 jq\, 
=0.025 and K c / = 0.^pM the relationship 

[Ca 2+ ] hasal ^ 07q2/qi implies that 



Pb = 



\D] 



[Ca 2+ ] bmal +K d 



[Ca + ]/, a5a /S0.014 J uAf which is small but not an unreasonable 
value. We proceed similarly for set (Hi). Namely, we use Eq. 12 
with y(l +qi + q 2 ) «8.18 on the points of Fig. 5C with (F}<22 
to get a range of possible values of — y 2 qiq 2 ^Ny. Using Eq. 8 to 
write y 2 q l q2(N) = y(F)q 2 /(pb + q2/qi) we find that pb is 
negligible with respect to q 2 /qi- This means that the basal 
fluorescence is dominated by that of the free dye molecules in this 
case. We then conclude from our estimates that, at basal 
conditions almost all of the fluorescence comes from the dye free 
molecules in the case of Rhod-2 while the contribution of the 
Ca 2+ -bound dye molecules is comparable to that of the free ones 
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in the case of Fluo-4. Assuming that qi/qi is of the order of the 
ratio of quantum efficiencies estimated for Rhod-2 in [24], 
<72/<7l =0.07, the relationship 

[Ca 2+ ] haml 



[CaD] b 



[Ca 



.2+1 



J hasal 



+ K d 



iqi/qi with K d ^2pM (the 



Pb-- 

value provided by Invitrogen) implies that 
[Ca + ] hasal < <0.15^Mwhich is consistent with the estimate 
obtained for the case of Fluo-4 ([C« 2+ ] Aaja ,s0.014jitM). Further- 
more, for the few points of set (Hi) for which we could estimate pb 
we obtained pb~0. 18^2/91 which implies [Ca 2+ \ hasal = 0.026^M, a 
value within the same order of magnitude as the one derived for 
Fluo-4. 

The second test consists of comparing the range of <A> that 
can be inferred from the experiments of Type I for cases (i) and 
(Hi) (Figs. 3A and 5A) with those that can be derived from the 
experiments of Type II (Figs. 3B and 5B). Particularly, from the 
experiments illustrated in Fig. 3B we can also estimate (N) in the 
observed region if we assume known values of p=Pb and qi/qi- 
The fit of the curve in Fig. 3B, which is not very good, estimates 
the prefactor of the nonlinear term, (Fy 2 , as 0.04. Thus, 

1 p + q\/q\ 

- — ; ^=0.04 according to Eq. 11 from which we 

< N >(p + q2M 2 

obtain <A>«140 if we use p b = 0.125 and <A>«250 if we use 
Pb = 0.03. These numbers are between 1.4 and 2 times larger than 
the largest one derived from Fig. 3A. The difference can be 
attributed to the low quality of the fit of Fig. 3B. Proceeding 
analogously for the case of set (Hi), we estimate <A r >= 115 from 
thefitofFig. 5B using/? =Pb = 0.05. This value is within the range 
of (Ny derived from Fig. 5A for the same case. The estimate 
becomes <A> = 78 for />/, =0.01 (which corresponds to 
[Ca 2+ ] hamt = 0.02/.iAf). This value is smaller than the lowest 
estimate derived from the fit of Fig. 5A and seems to indicate that 



[C, 



,2+1 



/ >0.02/.(M in the experiments probed in Fig. 5B. 



The third test consists of comparing the range of <(iV> values 
obtained for the three sets of concentrations among themselves. 
From the fit to Fig. 3A for the set of concentrations (i) we found 
41<<JV><65 using p/, = 0.125 (which corresponds to 
[Ca 2+ ] hasal = 0ApM using K d ^0.8fiM) and 113<<iV><177 
using pb = 0.03 (which corresponds to [Ca 2+ ] hasal = 0.02pM). The 
fit of Fig. 3B did not provide reliable estimate of <iV>. From the fit 
to Fig. 4A for the set of concentrations (it) we obtained 
39<<JV><98 using p=pb = 0.05 (which corresponds to 
[Ca 2+ ] hasal = 0ApM and K d ^2\iM) and 58<<iV><145 using 
/)/, = 0.01 (which corresponds to [Ca 2+ ] hasal = 0.02juM). These viuues 

are consistent with the ones derived for case (x) which have the same dye and EOTA 

concentrations. From the fit to Fig. 5A (Hi) we obtained 
98<<JV><177 for pb =0.05(which corresponds to 
[Ca 2+ ] hasal = QApM) and 145<<iV><260 for p h = 0M (which 
corresponds to [Ca 2+ \ hasal = 0.02pM). These values of <JV> are 
between 1.8 and 2.5 times larger in Fig. 5 A than in Fig. 4A which 
is consistent with the fact that the dye concentration in Fig. 5A is 
2.5 times larger than in Figs. 3A or 4A. 

Summarizing, our method passed successfully a set of self- 
consistency tests and allowed us to estimate y, q\ and qi at the 
standard illumination power (the latter assuming known values of 
qi/qi) an d a range of values of <iV> in the three examples probed 
in the present paper. 

Signal-to-noise ratios estimated with the method 

When we applied the method to estimate the signal-to-noise 
ratios of the experiments illustrated in Figs. 1, it again performed 



very well. More specifically, the method determined that the 
expected signal-to-noise ratio for the conditions of Figs. 1A,1C 
were similar while the one of Fig. 1 B was much smaller. This is 
consistent with the fact that we could not observe Ca 2+ puffs 
under the experimental conditions of Fig. IB, but we did observe 
them with the other two combinations. Thus, our method classifies 
the experiments of the example correctly in terms of signal 
detectability. Furthermore, the method can be used to determine 
to what extent changes in certain experimental parameters can 
lead to a noticeable improvement of the signal-to-noise ratio as 
illustrated in Fig. 8. This figure also shows how the method can be 
used to guarantee signal comparability between experiments 
performed with two different dyes when the dye used in a series 
of experiments needs to be changed because it is discontinued or 
replaced by an upgrade. 

Additional information that can be inferred with the 
method 

Having a realistic fluctuation model, on the other hand, allows 
us to go beyond a mere classification and draw additional 
information on the imaging experimental conditions. In particular, 
we can investigate in more detail what factors are most important 
for the differential ability of the different experimental settings to 
detect Ca 2+ elevations. We can do so in the example analyzed in 
this paper. More specifically, we can explore to what extent the 
different kinetics of the dyes and the different EGTA concentra- 
tions are responsible for the different types of behaviors observed. 
To explore the effect of the different dye kinetics on the observed 
images, we analyzed the blurred Ca 2+ -bound dye concentration 
obtained with the simulations. The obtained concentrations are 
shown in Fig. 6. A direct observation of this figure does not show 
significant differences in detectability among the three simulated 
situations that correspond to the set of concentrations (i) (A); (it) (B) 
and (Hi) (C). This coincides with a more quantitative comparison 
of the simulations. In particular, the ratio of the (blurred) Ca 2+ - 
bound dye concentration at the peak of the simulated signal over 
its basal value is approximately the same for all three figures. We 
investigated to what extent these results depended on the 
simplified model that we used for the simulations. In particular, 
we added a mobile buffer and explored a wide range of these 
parameters always obtaining the same behavior. The results of 
Fig. 6 and these additional studies imply that the differences in 
detectability observed between the cases (i) or (Hi) and those of 
case (ii) cannot be attributed to differences in the underlying 
Ca 2+ -bound dye distribution. It is only via the generation of noisy 
numerically simulated images by means of our quantified 
fluorescence fluctuation model that we can reproduce the 
detectability properties of the three experimental conditions 
probed in Fig. 1 . This also shows that the ratio given by Eq. 1 is 
not always a faithful reporter of the increment in the Ca 2+ bound 
dye concentration during signals with respect to the same 
concentration at basal conditions. 

The conclusion according to which the differences in detect- 
ability between the experiments of Figs, 1A and IB cannot be 
accounted for by differences in the underlyingCa 2+ -bound dye 
distribution led us to ponder the role of the Ca 2+ -free dye 
fluorescence on the images. To this end, we generated noisy 
images combining the fluctuation model and numerical simula- 
tions but without including the contribution of the Cfl 2+ -free dye 
molecules (i.e., setting q2 = 0). We show the results in Figs. 7D-F. 
Differences in detectability of the simulated signal are much less 
clear than in the case displayed in Figs. 7A-C. This observation is 
confirmed by a quantitative comparison. We show in the first 
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column of Table 5 the signal-to-noise ratios obtained when using 
q2 = 0. There we see that the ratio obtained for the case (it) is only 
slightly smaller than the one for case (i). Furthermore, the ratio for 
case (iii) is even larger than the one for case (i). This does not 
agree with what is observed in the Type 0 experiments performed 
for these sets (e.g., Fig. 1). We then conclude that noise by itself 
cannot account for the inability of Rhod-2 to report puff 
occurrences when used at low concentrations. The contribution 
of the Ca 2+ -free dye to the fluorescence is relevant for the 
differences in detectability observed when using Fluo-4 or Rhod-2. 
This is also apparent in the estimates of pi, when compared with 
<?2/<7l that can be derived combining the fits of Figs. 3A and 3C or 
5A and 5C. Namely, from Figs. 3A and 3C we obtained 
Pb = (0-l — 0.7)^2/?lfor the set of concentrations (i). This means 
that the contribution of the fluorescence of the Ca 2+ -bound dye 
can be comparable to that of the free dye under basal conditions 
for this setting. Combining the results of Figs. 5A and 5C we 
obtained Pb<qi/q\ (and pbO. IRqj/qi at most) for the set (iii). 
This means that basal fluorescence is mainly due to the free dye 
molecules for this experimental setting. 

The different values of [EGTA] used in the experiments of 
Figs. IB and 1C also play a relevant role on the better detectability 
properties of set (iii) when compared with those of set (it). This is 
reflected in the different behaviors that the <r F vs(Fy curves 
displayed for the two experimental settings when changing the 
illumination power (Figs. 4B and 5B). Namely, while this curve has 
a parabolic shape in Fig. 5B, it remains a straight line in Fig. 4B. 
This means that the nonlinear part of the a F vs <.F> curve is not 
reached for the laser powers probed in the case (it). According to 
Eq. 1 1 the nonlinearity of the oj-vs (F") should become noticeable 

(Fy 2 p^rQ 2 /^ 2 

when a = y(F}<b= — — ; — — Kr. Inserting Eq. 8 into this 

inequality we obtain q\ >c = ^—5^1. The quantity c>\ since 
P + qi/qi 

<72/<7l<E Furthermore, it is c=l if p>>qijq\ while it is 
c = q\fq2>>\ for p< <qijq\. This implies that the smaller 
/>with respect to qi/qi, the larger q\ needs to be for the nonlinear 
term of Eq. 1 1 to be noticeable. Therefore, the different behavior 
of Figs. 4B and 5B could be attributed to a smaller value atp=pb 
in the experiments of Fig. 4B. This, in turn, could be due to the 
larger amount of EGTA used in those experiments compared to 
the ones of Fig. 5B. This analysis together with the fact that no 
increment in fluorescence was observed upon Ca 2 + microinjec- 
tion for the set of concentrations (it) is an indication that EGTA is 
capturing Ca 2+ more efficiendy than the dye for this experimental 
setting. This explains why a larger [Dye] and smaller [EGTA] is 
necessary to report the occurrence of localized signals using the 
dye, Rhod-2, with similar detectability levels as when using Fluo-4. 
This analysis not only illustrates the possible applications of our 
method. It also indicates that, in certain cases, it is not necessary to 
compute the signal-to-noise ratio by means of numerical 
simulations to determine the inadequacy of certain combinations 
of dye and EGTA concentrations to observe localized 
Ca 2 + signals. In particular, it seems that we should not expect to 
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Final remarks 
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of Ca 2+ -bound and the Ca 2+ -free dye molecules, FLIM would 
give an accurate in situ measurement of pb and of qi/qi- Our 
method would then provide estimates ofy, <JV>and q\ and, by 
means of numerical simulations, signal-to-noise ratios and 
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